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ABSTRACT 

We combine A^-body simulations of structure growth with physical modelling of galaxy 
evolution to investigate whether the shift in cosmological parameters between the 
1-year and 3-year results from the Wilkinson Microwave Anisotropy Probe affects 
predictions for the galaxy population. Structure formation is significantly delayed in 
the WMAP3 cosmology, because the initial matter fluctuation amplitude is lower on 
the relevant scales. The decrease in dark matter clustering strength is, however, almost 
entirely offset by an increase in halo bias, so predictions for galaxy clustering are barely 
altered. In both cosmologies several combinations of physical parameters can reproduce 
observed, low-redshift galaxy properties; the star formation, supernova feedback and 
AGN feedback efficiencies can be played off against each other to give similar results. 
Models which fit observed luminosity functions predict projected 2-point correlation 
functions which scatter by about 10-20 per cent on large scale and by larger factors 
on small scale, depending both on cosmology and on details of galaxy formation. 
Measurements of the pairwise velocity distribution prefer the WMAP1 cosmology, but 
careful treatment of the systematics is needed. Given current modelling uncertainties, 
it is not easy to distinguish the WMAP1 and WMAP3 cosmologies on the basis of 
low-redshift galaxy properties. Model predictions diverge more dramatically at high 
redshift. Better observational data at z > 2 will better constrain galaxy formation and 
perhaps also cosmological parameters. 
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1 INTRODUCTION 

Our current understanding of the cosmic evolution is 
based heavily on measurements of the Cosmic Mi- 
crowave Background (C MB). The discovery of the CMB 
|Penzias fc Wilson! Il96sh led to general acceptance of the 
Hot Big Bang theory, and further support came from the 
subsequen t demonstration of its near-perfect black-body 
spectrum (Ma ther et al.l ll990). The detection and continu- 
ing refinement of measures of angular structure in the CMB 
confirmed theoretical predictions for the growth of struc- 
ture in flat cosmolog ies dominated by non-baryonic Cold 
Dark Matt er (CDM) (|Smoot et al.lll990l : lde Bernardis et al.1 
l200d . l2002h . The power spectrum of this structure encodes 
information about the values of the cosmological parame- 
ters, although degeneracies prevent a precise determination 
of all parameters from CMB data alone. By including con- 
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straints from other kinds of data, it becomes possible to 
constrain many cosmological parameters quite tightly. Dur- 
ing the last decade, a growing body of such measurements 
have ushered in a new era of 'precision cosmology'. 

The first-year data fro m the Wilkinson Microwave 
Anisotropy Probe (WMAP) (|Spergel et al.ll2003f ) did much 
to establish ACDM, a flat CDM model with a cosmolog- 
ical constant, as the standard model for structure forma- 
tion. In addition, they provided apparently precise esti- 
mates for a number of cosmological parameters. Two fur- 
ther years of WMAP data have significantly refined these 
estimat es, leading to notic eable shifts in some of the best 
values (jSpergel et alj|2007f ). The most important differences 
are a reduction in optical depth to the last scattering surface 
(r), a lower value for the amplitude of matter fluctuations 
on 8/i _1 Mpc scale (erg), a reduction of the scalar spectral 
index on primordial perturbations (n), and a lowering of the 
total matter density (S7 m ). 

Several studies prior to the latest WMAP results 
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suggested a lower value for as- These loo ked at evolu- 
tion i n the abundance of g alaxy clusters (|Borgani et al.l 
l200ll : ISchuecker et al.l 12003m . and at constraints from 
the abundance and clustering of low-redshift galaxies in 
combination with observed cluster mass-to -light ratios 
l|van den Bosch et all 120031 : iTinker et all 120051 ) . Parameter 
estimates from these methods are partially degenerate in 
as and Sl m and other studies, notably of cosmic shear, 
have t ended to give conflicting indications (e.g.lMassev et al.l 
120071 ; iBeniamin et alj 120071 ; Irloekstra et alj"l2006l) . Recent 
studies on giant arc statistics (|Li et aL 20061 ) also suggest 
that low values of as may be difficult to reconcile with the 
observations. It is only the shrinking of the allowed param- 
eter region forced by the new WMAP data that has per- 
suaded much of the astronomical community to prefer a 
"standard" model with a lower value of as- It should be re- 
membered, however, that including parameters beyond the 
usual minimal set significantly relaxes constraints, so that 
both first- and third-year WMAP parameter sets should 
probably be treated as plausible. 

Modifications of the cosmological parameters of the 
kind discussed above can significantly affect the formation 
and evolution of structure. Recent studies have discussed the 
implications of WMAP third-year results for the formation 
of the first stars and black holes, and for the reionization 



Table 1. Cosmological parameters of the two simulations used 
in this study. fi ra , Q^, f2j, are the density of matter, dark energy 
and baryons respectively, as and n are the amplitude of the mass 
density fluctuations, and slope of initial power spectrum. The 
Hubble constant is parameterised as Hq = 100 h km s~ 1 Mpc _1 . 



of th e intergalact i c medium llGao et alj 20061; Alvarez et al.l 
120061 ; |Popal l2006l ; iLewis et al.ll2006l ; llliev et alj|2007h . rjncT 
ing significantly later formation and reionization than pre- 
viously thought. 

The present paper investigates the impact of this de- 
lay on the observed properties of galaxies, both at low and 
at high redshift. We combine high-resolution iV-body simu- 
lations with semi-analytic modelling te chniques to simulate 
the evolution of the galaxy population llSpringel et al.ll2005l ; 
ICroton et al.ll2006l ; lDe Lucia fc Blaizotj|2007h . Our paper is 
organised as follows: the iV-body simulations and the semi- 
analytic modelling assumptions are described in sections 
and [3] In section [4] we discuss the formation of dark mat- 
ter structures in cosmologies with first-year and third-year 
WMAP parameters, while in section [S] we compare simu- 
lations of the evolution of the galaxy population in these 
two cosmologies. Finally, in section[S]we summarise and dis- 
cuss our findings. For the convenience of the reader we also 
provide an appendix containing an analytic comparison of 
evolution in the abundance and clustering of dark halos in 
our tw o cosmologies, based on the formulae of lMo fe White! 
I2002T ) 



2 THE SIMULATIONS 

We have carried out two simulations of the growth of struc- 
ture in a ACDM Universe using the cosmological parameter 
sets listed in Tab. [T] The cosmological parameters used for 
our WMAP1 simulation are derived from a combination of 
first-year WMAP results dSpergel et al.ll2003ri wit h the 2dF- 
GRS galaxy power spectrum IjColless et al J 120011 ). and cor- 
respond to those us ed for the Millennium Simulation (MS; 
ISpringel et aLl feoOS). Our WMAP3 simulation adopts cos- 
mological parameters derived from a combination of third- 
year WMAP data on large scale, and Cosmic Background 
Imager (CBI) and extended Very Small Array (VSA) data 





WMAP1 


WMAP3 




0.25 


0.226 


n A 


0.75 


0.774 


n b 


0.045 


0.04 




0.9 


0.722 


h 


0.73 


0.743 


n 


1 


0.947 



on small scale l|Spergel et al.l [2007l "l . We note that, among 
the different sets of cosmological parameters consistent with 
the third-year WMAP data and other observational data, 
we have chosen to look at the one which differs the most 
from the parameters of our WMAP1 simulation. 

The most significant difference between the two sets 
of parameters listed in Table [T] is in the lower value of 
a 8 adopted for the WMAP3 simulation. Our WMAP1 and 
WMAP3 cosmologies also differ in the scalar spectral index 
of primordial density perturbations (n = 1 — > 0.947), and in 
the matter density (f2 m = 0.25 — » 0.226). As noted above, 
a number of recent studies have shown that this change in 
cosmology results in a significant delay of structure forma- 
tion in the WMAP3 case in comparison to WMAP1. We will 
discuss this in more detail in the following section. 

The numerical parameters used for our simulations are 
listed in Table [2] The mass and force resolution are the same 
as used for the MS, while the volume is a factor 64 smaller. 
We have stored the data at the same 64 output times as in 
the MS. These are approximately logarithmically spaced be- 
tween z = 20 and z — 1 and approximately linearly spaced in 
time thereafter. Friends-of-friends (FOF) group catalogues 
were computed on the fly for each sna pshot, and the al- 
gorithm SUBFIND l|Springel et all bOOll l was employed to 
decompose each FOF group into a set of disjoint substruc- 
tures. As in lSpringel et all l|2005h . only substructures which 
retain at least 20 bound particles after a gravitational un- 
binding procedure are considered to be genuine substruc- 
tures. Substructure catalogues are then used to construct 
detailed merger history trees that provide the basic input 
needed for the semi -analytic mode l desc ribed in the next 
section. We refer to ISpringel et all (|2005T ) for more details 
of the merger-tree construction. Both our simulati ons were 
run u sing the tree-based parallel code GADGET2 (|Springell 
2005). The initial power spectra w ere generated using CMB- 
FAST (|Seliak fc Zaldarria ga 1996) with the cosmological pa- 
rameters listed in Tab. 1 as input. The Fourier modes of the 
initial density field in the two cases were identical except for 
the amplitude adjustment needed to reproduce the correct 
power spectrum. Thus structures correspond closely in the 
two cases. 

Fig. [1] shows the power per decade in the linear initial 
conditions of our WMAP1 (blue) and WMAP3 (red) simu- 
lations. The significantly lower value of as in the WMAP3 
case reduces the overall normalisation. In addition, the red 
tilt in the primordial power spectrum index makes the differ- 
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Table 2. Numerical parameters of our two simulations. L is the 
box size, n p is the particle number, e is the softening length, and 
m p is the particle mass. 



L n p e m p (WMAPl) m p (WMAP3) 

125MPC/1" 1 540 3 5kpc 8.61e8M h" 1 7.78e8M h' 1 



0.10 1.00 
k[h/Mpc] 



Figure 1. The power spectra of our two simulations: WMAP1 
(blue solid) and WMAP3(red dashed). These were generated us- 
ing CMBFAST with the cosmological parameters listed in Tab.[T] 
The quantity plotted is A 2 (fc) oc k 3 P(k). 



ence greater on small scales than on large scales. Indeed, on 
the scales responsible for the low £ microwave background 
fluctuations the two power spectra are almost the same, re- 
flecting the fact that they are both required to fit observed 
CMB fluctuations on these scales. Conversely, the suppres- 
sion of power is particularly significant on the small scales 
responsible for the formation of the first nonlinear struc- 
tures. 



3 THE SEMI-ANALYTIC MODEL 

In this p aper we use the galaxy formation model de- 
scribed in lDe Lucia fc BlaizotJ (|2007h which is based on the 
WMAP1 cosmology and builds o n previous work by the 
"Munich" galaxy formation group (Kauffma nn et al. 1999a 



Sprmgelet all 2001 



2005; ICroton et al.l 



De Lucia et all 12004 ISpringel et al l 
Although not in perfect agree- 



2006) 



ment wit h all aspects of the local galaxy population (see, for 
example. IWeinmann et al.l (120061 1). this model does quite a 
good job of reproducing the observed relations between ste l- 
lar mass, gas mass, and metallicity |De Lucia et al.l I2004T ) . 
the o bserved luminosity, colour, and morphol ogy distribu- 
tions (|Croton et al.ll2006l ; |Pe Lucia et all 120061) and the ob - 
served two-point correlation functions (|Springel et al.l l2005). 



licit zbichler fc White! <|2007h have recently shown that it also 
agrees reasonably well with the observed galaxy luminosity 
and mass functions at high redshift. We refer the reader to 
the original papers for a full description of the numerical im- 
plementation, and of the physical processes modelled. In the 
following, we summarise briefly the treatment of those phys- 
ical processes for which we needed to change the efficiency 
parameters in order to maintain agreement with observa- 
tions of the local galaxy population when we switch from 
WMAP1 to WMAP3 (see Tableland Sec. EJ- 

In the semi-analytic model we use in this work, star 
formation is assumed to occur at a rate given by: 



TO* = QSF(m co l d — m cri t)/tdyn,di S 



(1) 



where m co id and tdyn.disc are the cold gas mass and the dy- 
namical time of the galaxy, defined as the ratio between 
the disk radius and the virial velocity, m cr i t c orresponds 
to a critical value for t he gas surface density ( Kennicuttl 
1 19981 : lKauffmannlll996l ; iMo et al.lll998D . and qsf controls 
the efficiency of the transformation of cold gas into stars 
when the gas surfac e density is above the critical value. (See 
ICroton et al.l (|2006h for more detailed descriptions of the im- 
plementation of this process and of the feedback processes 
described below.) 

Massive stars explode as supernovae shortly after a star 
formation event and are assumed to reheat a gas mass that 
is proportional to the mass of stars formed (based on the 
observations of iMartinll 19991 ): 

Am rc hcatod = EdiskAm,, (2) 

Again following ICroton et all (|2006l ). we write the energy 
released by an event which forms a mass Am* in stars as: 



A£ SN = 0.5e ha loATO.Vs 



(3) 



where 0.5V S 2 N is the mean supernova energy injected per 
unit mass of newly formed stars, and ehaio represents the 
efficiency with which this energy is able to convert cold in- 
terstellar medium into hot, diffuse halo gas. The amount 
of gas that leaves the dark matter halo in a "super-wind" 
is determined by computing whether excess SN energy is 
available to drive the flow after reheating of material to the 

halo virial temperature. 

As in iKauffmann fc HaehnelH l|2000h . black holes are 
formed and fuelled during mergers: 



Am B H 



/bH™c 



1 + (280kms- 1 /Kir) 2 



(4) 



We assume here that black holes grow during both major 
and minor mergers, and that the efficiency of gas accretion 
onto the black hole scales with the baryonic mass ratio of 
the merging galaxies: 



/bh = /bh x (m sat /m ccntra i) 



(5) 



This is the primary process driving the growth of the total 
mass in supermassive black holes. Individual black holes can 
also gain mass through merging when their host galaxies 
merge. 



Finally, we use the model bv lCroton et al. (2006) to de- 
scribe heating by centrally located AGN in massive groups 
and clusters. This process is assumed to be associated with 
"radio mode" outflows which suppress cooling flows and thus 
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the condensation of gas onto the central galaxies. The pro- 
cess is assumed to occur whenever a massive black hole finds 
itself at the centre of a static hot gas halo. During this phase, 
the accretion rate onto the central supermassive black hole 
is taken to be: 

( ™BH \ //hot\ / Kir \ 3 ,„v 

mBH ' R = KAGN \WM7, ) {^IJ UoOkms-0 ' (6) 

where mBH is the black hole mass, /hot is the fraction of the 
total halo mass in the form of hot gas, Kir is the virial veloc- 
ity of the halo, and kagn is efficiency parameter with units 
of Af©yr . The energy released durin g this accretion pro - 
cess is used to reduce the cooling flow. ICroton et alj (|2006l ) 
showed that this results in complete suppression of cooling 
in relatively massive haloes and groups. The process starts 
being effective at a mass scale that evolves as a function of 
redshift as shown in their Fig. 7. 

In our simulation scheme, haloes (and the galaxies 
within them) are followed even after they are accreted onto 
larger systems. The dynamics of such a satellite subhalo 
is followed explicitly by the N-body simulation until tidal 
stripping causes i ts mass to fall below the resolution limit 
of the simulatio n (|Ghigna et alj|200d : I5e Lucia et alj|2004t 
iGao et ahl Eool). When this happens, we estimate a sur- 
vival time (tmerge) for the associated galaxy using its cur- 
rent orbit and the classic al dynamical friction formula of 
iBinnev fc Tremaind (|l987j ). Once this time has elapsed, the 
galaxy is assumed to merge onto the central galaxy of its 
current halo. While it still survives it is assumed to follow 
the particle which was the most bound particle of the sub- 
halo a t the last time it was identified. IDe Lucia fc Blaizotl 
|2007l ) found that increasing the merging time by a factor of 
2 slightly improves the fit to observed luminosity function. 
Such a n increase has other effects which Il5e Lucia fc Blaizotl 
120071 ) did not study, for example, it increases the ampli- 
tude of small-scale galaxy correlations by about a factor of 
2 at r < 100 kpc. Others authors have claimed that this 
merging time should be effectively set to zero, so that such 
"orphan" galaxi es lose their identit y at the same time as 
their subhaloes l|Conrov et al.l 120071 ). Here we leave this is- 
sue for detailed study in future work and simply consider 
the pre-factor to be used for tm CVgc as a free parameter. 



4 FORMATION OF DARK MATTER 
STRUCTURES 

Figure [5] shows the dark matter distribution at z = within 
a slice 20 Mpc h~ x thick cut from the full volume of our sim- 
ulations. The WMAP1 model is on the left and the WMAP3 
model on the right. The projection is colour-coded by den- 
sity and provides a visual illustration of the delay of struc- 
ture formation in the lower erg model. Although the overall 
structure is very similar, it is clear that massive haloes lying 
at the nodes of the cosmic web are in a more advanced stage 
of merging in the WMAP1 case than in the WMAP3 case. 

The differing fluctuation amplitudes of the two simu- 
lations translate into a systematic difference between their 
halo mass functions. This is illustrated in Fig. [3] where we 
compare the cumulative number density of halos for the two 
cases (blue for WMAP1, red for WMAP3). In addition, we 
show the corresponding function for the much larger Mil- 



lennium Simulation (green) which used our WMAP1 cos- 
mological parameters. It agrees very well with the smaller 
WMAP1 simulation of this paper. At z = 0, the most mas- 
sive haloes in our WMAP3 simulation are about 1.5 times 
less massive than their counterparts in the WMAP1 simula- 
tion. The number of haloes more massive than 10 13 M Q h' 1 
is ~ 25 per cent smaller in the WMAP3 case than in the 
WMAP1 case. These differences increase at higher redshift. 
At z ~ 6 (i.e. at the end of the reionization epoch) the most 
massive haloes in the two simulations differ by about a factor 
of 5 in mass. 

We recall that the main differences between the cosmo- 
logical parameters used in the two simulations are the lower 
value of erg and the redder primordial power spectrum index 
n in the WMAP3 case. As shown in Fig. [T] these combine 
to produce a substantial suppression of small-scale power. 
This has important implications for the formation of the 
first objects and for the star formation history at early times. 
IGao et all (|2006l ) studied this problem using high-resolution 
cosmological simulations with a smoothed particle hydrody- 
namics treatment of baryonic processes. They found a much 
lower abundance of potentially star-forming halos at high 
redshift for WMAP3 than for WMAP1. This reduction in 
the number of 'm ini-haloes' at z ~ 10 was also discussed by 
iReed et al.l (|2007l ) using analytic models to explore the de- 
pendence on cosmological parameters. For fluctuation am- 
plitudes at the WMAP1 level, very efficient production of 
UV radiation is needed to reionize the intergalactic medium 
by z ~ 15, as required by the WMAP1 value for th e electron 
scattering optical depth r fe.g. lCiardi et al1l2003l) . Interest- 
ingly, the delay in structure formation in a WMAP3 Uni- 
verse is such that equally efficient UV production is needed 
to reionize by z ~ 10, as inferred from the WMAP3 value of 
r (|Alvarez et al.ll2006l ). 

Fig. [4] shows z = autocorrelation functions in our 
WMAP1 and WMAP3 simulations for the dark matter (in 
the top panel) and for all dark matter substructures iden- 
tified by SUBFIND (in the bottom panel; this corresponds 
to all self-bound (sub)haloes with more than 20 simulation 
particles). The difference in fluctuation amplitude causes the 
correlation strength to be almost a factor of two smaller in 
the WMAP3 case than in the WMAP1 case on scales below 
a few Mpc. Curiously, however, this decrease is almost com- 
pletely compensated by an increase in the (sub)halo bias, 
so that the 2-point correlation functions for subhaloes are 
almost identical in the two cases. There is a slight residual 
offset on scales 0.5Mpcrt _1 < r < 5Mpc/i _1 . As we will see 
below, this results in very similar galaxy correlation func- 
tions being predicted in the two cases. 



5 GALAXY FORMATION 

In this section we analyse galaxy populations simulated us- 
ing the semi-analytic framework presented in Sec. [3] Table[3] 
lists the combinations of semi-analytic parameters for which 
we will show detailed results. Mod el A is the parameter set 
used bv lDe Lucia fc Blaizotl ([2007) for the Millennium Sim- 
ulation. As expected, this also gives good agreement with 
observation for our WMAP1 simulation. Models B and C 
are parameter combinations that, as we show below, pro- 
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Figure 2. Rcdshift zero distribution of dark matter within a slice of thickness 20Mpch through our two simulations, WMAP1 (left) 
and WMAP3 (right). 
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Figure 3. Cumulative halo number density as a function of halo 
mass and redshift. Results for the WMAP3 and WMAP1 simula- 
tions are represented by red and blue curves respectively. Results 
for the Millennium Simulation are shown in green. The halo abun- 
dance is multiplied by M 2 in order to reduce the dynamic range 
of the ordinate. 
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duce a similar agreement with local data for the WMAP3 
cosmology. 

Model B has the same star formation efficiency as model 
A («sf), but lower supernova feedback efficiency (ehaio and 
tdisk), and lower AGN feedback efficiency (kagn)- For this 
model, we have also eliminated the pre-factor of 2 which 



Figure 4. Autocorrelation functions at z = for dark matter 
(top panel) and for all subhaloes with at least 20 particles (bottom 
panel). Blue solid and red dashed curves show results for our 
WMAP1 and WMAP3 simulations respectively. 
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Table 3. Principal parameters of our galaxy formation models 
(see text for details). 



ABC 



asF 


0.03 


0.03 


0.07 


Ehalo 


0.35 


0.2 


0.28 


edisk 


3.5 


3 


4.5 


KAGN 


7.5e-6 


6.5e-6 


1.2e-5 


/bh 


0.03 


0.03 


0.05 


Emerge 


2 X ^friction 


^friction 


2 X tfriction 


Best for 


WMAPl 


WMAP3 


WMAP3 



|Pe Lucia fc BlaizotJ l|2007l ') introduced in the definition of 
the merging time. In contrast, in model C, we double the 
star formation efficiency (relative to model A) in order to 
compensate for the delay in structure formation, but this 
must be compensated by much higher feedback efficiencies 
(both from supernovae and from AGN) to prevent the over- 
production of stars at late times. For this model we keep 
the pre-factor of 2 in the definition of the merging times of 
satellite galaxies. We also increase the efficiency of accre- 
tion onto black holes during mergers (which increases the 
effectiveness of the 'radio' mode - see Eq. 6). In the follow- 
ing, we will show that models B and C give similar z = 
results for the WMAP3 cosmology, and that these resem- 
ble the results of model A for the WMAPl simulation. This 
shows that there are at least two independent way to 'com- 
pensate' for the delay in the structure formation within this 
physical framework: we can decrease feedback on all scales 
while keeping the same star formation efficiency (model B), 
or we can increase both the star formation efficiency and the 
feedback efficiency on all scales (model C). 

In the following sections we discuss how the predictions 
of our three models compare to observational data both in 
the local Universe and at high redshift. 

5.1 Low-redshift Luminosity Functions 

Figure [5] compares observational estimates of the K-band lu- 
minosity function of nearby galaxies to predictions from our 
three galaxy formation models (differentiated by colour) ap- 
plied to each of our two simulations (differentiated by line 
type) . In both cosmologies models B and C give almost iden- 
tical results, while model A predicts fewer galaxies at bright 
luminosities (Mk < —22) by a factor that varies from 1.3 to 
3. The inset in Fig. [S] repeats the figure showing only results 
from model A applied to our WMAPl simulation and from 
models B and C applied to our WMAP3 simulation. This 
shows all three models to produce similarly good fits to the 
observations around and above the 'knee'. All three models 
show an excess of galaxies fainter th an Mk — 51og/i ~ — 22 
with respect to the lCole et al.l (|200ll ) d ata, although t hey li e 
below the luminosity function given bv lHuang etall l|2003h . 
This possible excess is more pronounced in model A. In the 
following we will only discuss results from our 'best' mod- 
els (i.e. model A for our WMAPl simulation and for the 
Millennium Simulation, models B and C for our WMAP3 
simulation; in the following we denote these as Al, MSA, 
B3 and C3). 

Figure[6]compares model galaxy luminosity functions in 




-26 -24 -22 -20 -18 

M k - 5loc 10 h 



Figure 5. K-band luminosity function for the three galaxy for- 
mation models and the two cosmologies simulated in this study. 
Results from models A, B, and C arc shown in blue, red and green 
respectively. Solid and dashed lines correspond to our WMAP3 
and WMAPl simulations respectively. The black sy mbols with 
error bars show the observational determination by ICole et al.l 
(l200lh. while the black dotted curve shows the measurements of 
iHuane et al. | d2003h . The inset repeats the figure but shows re- 
sults only from model A applied to WMAPl and from models B 
and C applied to WMAP3. 



1 cr 1 




-24 -23 -22 -21 -20 -19 -18 -17 

M .1,r 

Figure 6. Galaxy luminosity function at SDSS ro.l band. 
Coloured lines show predictions from the three "good" models 
of this paper. We also plot results from model A applied to the 
Millennium Simulation. The black curve with an e rror band cor- 
respon ds to the observational measurement from iBlanton et al.l 
J2003h . 
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the Sloan Digital Sky Survey (SDSS) r-band to observational 
measurements from the SDSS itself. Here the observational 
uncertainties are much smaller than for the K-band lumi- 
nosity function of Fig. [5] Since the obs erved absolute mag- 
nitudes quoted bv lBlanton et all ((2003) are band-shifted to 
2 = 0.1, we also correct our simulated r-band absolute mag- 
nitude to this redshift using their K-correction code. Fig. [6] 
shows all our "good" models to agree very well with the ob- 
servational data, particularly around the knee of the lumi- 
nosity function ( — 22 < M r < —19). All models overpredict 
the abundance of faint galaxies, although the problem ap- 
pears significantly less dramatic here than in Fig. [5] In all 
models the rarest and most luminous galaxies are also too 
bright by 0.2 to 0.6 magnitudes. This problem may be ,at 
least partially, due to our assigning light to these objects 
which should be part of the intracluster light of the groups 
or clusters of which they are the central galaxies (see also 
IConrov et all l|2007l )). 

5.2 The Tully-Fisher Relation 

In Fig. [7] we compare the Tully-Fisher relation for our 
model galaxies to t he observational determination by 
IGiovanelli et~ai1 < |l997h . Green symbols show all model cen- 
tral galaxies with 1.5 < Ms^ulge — Me.total < 2.6 where 
Mb, bulge and Ms.totai are the absolute magnitude of the 
bulge and the total B-band magnitude, respectively. This 
selection is made in order to isolate Sb/ c spirals as in the 
observed sample of lGiovanelli et al.l (|l997l ) . The mean obser- 
vational relation is shown by a solid blue line in the figure, 
with the corresponding la scatter indicated by the dashed 
lines. The relation of Giovanelli et al. is already corrected for 
internal extinction. We therefore do not include dust effects 
when predicting I-band magnitudes for this plot. Red lines 
in each panel show a linear fit to the model results. The top 
panel of Fig.[7]demonstrates that model Al reproduces both 
the slope and the zero-point of the observed relation (as also 
shown in lCroton et alj|2006h . Models B3 and C3 exhibit a 
substantially brighter zero-point than observed. This occurs 
because fitting the observed luminosity function within the 
less evolved halo mass function of the WMAP3 cosmology 
(see Fig. 2) requires us to put galaxies of given luminosity 
at the centre of lower mass halos than in the WMAP1 case. 

Finding theoretical models which simultaneously fit 
both the observed luminosity functions and the ob- 
served Tully-Fis her relation has always been difficult (e.g. 
ICole et alJl200Ch . The results in Fig. [TJ show that this re- 
mains true. Note, however, that it is uncertain how best to 
extract quantities from our models for comparison with the 
observational data. For example, we assume here that W 
(the measured HI linewidth) is twice the maximum circular 
velocity of the galaxy's dark halo, as measured directly from 
the simulation and without any correction for the effects of 
the baryonic mass of the galaxy. This is a rough assumption 
which may be systematically in error. In particular, if galaxy 
formation results in a maximum disk rotation velocity which 
is significantly and systematically above this value, the pop- 
ulation of galaxies in our WMAP3 model could be recon- 
ciled with the observations. Changing the transformation 
between these two velocities can change the slope and zero- 
point, or introduce curvature into the model Tully-Fisher 
relation. Semi-analytic models provide robust estimates for 
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Figure 7. Tully-Fisher relation for spiral galaxies in our three 
models. We show only central galaxies with 1.5 < Mg bulge — 
Mg total < 2-6, and we approximate the observed HI linewidth by 
twice V m ax for the dark halo, as measured directly in our simula- 
tions. Red lines are linear fits to the model result s. The blue solid 
line s hows the mean observational relation by IGiovanelli et all 
lll997l) , with blue dashed lines indicating the scatter in this rela- 
tion. 



the locations, velocities and global properties (luminosity, 
colour, stellar mass, etc,) of galaxies, but become less reli- 
able when quantities (such as W) which depend sensitively 
on the internal structure of the galaxies must be pred icted 
(see also the discussion in ISomerville fc Prim ack 1999). 

5.3 Mass-to-Light Ratios of Clusters 

In Fig. [8] we plot SDSS ro.i band mass-to-light ratios 
(M /L rt o.i) for individual clusters with M ^ 1O 14 M0 (sym- 
bols) and the running median value over the full resolved 
mass range (solid lines). Both M and L r ,o.i are computed 
within -R200, defined as the radius within which the mean 
mass overdensity is 200 times the critical value. Blue, red 
and green are used for models Al, B3, and C3 respectively. 
Cyan is used for the model MSA which provides a larger 
number of clusters than the smaller simulations used in our 
study. The black horizontal line and hatche d area show the 
region occupied by the observational data in lCarlberg et all 
(1996, see Tinker et al. 2005 for details on the conversion to 
the SDSS ro.i band). Our models exhibit a very weak mass 
dependence over the observed mass range and agree reason- 
ably well with the observational measurements. The differ- 
ences between the two cosmological models, and between the 
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two galaxy formation models for same WMAP3 cosmology 
are sm all, even at small masses where Ivan den Bosch et all 
(|2007l ) found a more pronounced decrease in the mass-to- 
light ratio for WMAP3 cosmology with respect to WMAP1. 

It is interesting to compare our results with those 
based on halo occupatio n distribution (HOD) models. 
Ivan d en Bosch et al . (2003) found that their conditional lu- 
minosity function (CLF) model was unable to fit the ob- 
served mass-to-light ratios of clusters in a cosmology with 
WMAP1 parameters, and they argued strongly in favour 
of a cosmology with a lower value of as, similar to that 
subsequently preferred by WMAP3. Once their models are 
constrained to fit the observed luminosity and correlation 
functions of galaxies, t hey find (M / L) to vary roughly 
as erf on cluster scales. iTinker et al.l l|2005l ) found a sim- 
ilarly strong dependence of the mass-to-light ratio on <rg 
using an HOD model which differed in detail and which 
they constrained with different observational data. M/L 
predictions from these models are shown in Fig. [8] by 
the horizontal dashed and dotted lines (blue and red are 
used for WMAP1 and WMAP3 respectively). Predictions 
from Ivan den Bosch et all (2003) were obtained by con- 
verting their B-band estimates to the SPSS rp.i band us - 
ing the conversion factor adopted in ITinker et al] (|2005h . 
These authors convert the mean B-band mass-to-light ra- 
tio < M/L B >= 363/i(M/L) of ICarlberg et all (|l996l ) to 
< M/Lrfl.i >= 359/i(M/ L) fl in the SPSS r .i band Thus 
we multiply the results of Ivan den Bosch et al] |2003) by a 
factor of 359/363 = 0.988. Fig. [8] shows that the mass-to- 
light ratio dependence on <rg in our models is much weaker 
than in the two HOP models mentioned above. This is sur- 
prising since our models are also a good fit to the observed 
galaxy luminosity functions, and fit observed correlation 
functions moderately well, at least on large scale (see be- 
low). 

In Fig. [8] we also show as a black solid line the average 
mass-to-light ratio predicted by the ph ysically based HOP 
model presented in IWang et al] (J2006) . These authors use 
the same positions and velocities for galaxies as our semi- 
analytic model (taken directly from the Millennium Simu- 
lation) but adopt simple parametrized functions to relate 
the star formation rates of galaxies to the mass and location 
(satellite or central) of their parent subhalo. The parameters 
of these functions are then minimized in order to optimize 
the fit to the observed luminosity function and luminosity- 
dependent correlation functions. The cluster mass-to-light 
ratios predicted by this model are higher than those in any 
of the semi-analytic models in our study, and are r eason - 
abl y close to the predic tions of Ivan den Bosch et all (|2003l ) 
and ITinker et al] 1120051 ) for a WMAP1 cosmology (the cos- 
mology adopted in the Mille n nium Simulation). Compar- 
ing the model of I Wang et~ail l|2006h in detail to the MSA 
model, we find that the larger mass-to-light ratio on clus- 
ter scales is due to systematic lower central galaxy lumi- 
nosi ties (by a factor ~ 2) in the HOP model. We note 
that lPe Lucia fc Blaizotl (|2007T ) found reasonable agreement 
when comparing this same MSA model to the observed 
K-band magnitudes of brightest cluster galaxies (BCGs). 
We note also that precise measurements of BCGs luminosi- 
ties are complicated by the intrinsic diff iculty in separat- 
ing the contribution of the cP envelope l|Schombertl 1 19881 ; 
iGonzalez etail 120051 ; IZibetti et al1l2005l ). The three HOP 



1 B88 

800 
700 
600 
© 500 



- 300 

d 
O 

-f 200 





















/ yj 




/ //■ — w jSitefi-i 

- 7-jrf J^r- ~ W7 1 




y A1 




B3 




C3 




MSA 




MSA-HOD 



1Q 12 1Q 13 1Q 14 1Q 15 

M 200 [M ] 

Figure 8. Mass-to-light ratios in the SDSS ro.i band as a func- 
tion of halo mass. Symbols show results for the 48 (20) clusters 
with mass M > 1O 14 M in our WMAP1 (WMAP3) simulation. 
Solid lines show a running median over a wider range of halo 
masses. Blue, red, green and cyan are used for models Al, B3, C3 
and MSA respe ctively. The black solid curve refers to the HOD 
model of I Wang et al.l ((200G) applied to the Millennium Simulation 
(MSA-HOD). The black dashed line and the hat ched area show 
the re gion occupied by the observational data of ICarlberg et aTl 
(1996). The horizontal dotted and dashed li nes show predictions 
for M /L of clusters M >10 14 M Q from | van den Bosch et all 
(2003) and Tin ker et al] (|2005| ) respectively. In each case blue 
refers to the WMAP1 and red to the WMAP3 prediction. 

models shown in Fig. [5] exhibit a stronger dependence on 
cosmology (in particular as), but also a larger dependence 
on modelling details than do our semi-analytic models. It 
seems, therefore, that a good understanding of the con- 
sequences of modelling assumptions is needed before firm 
conclusions can be drawn about as or other cosmological 
parameters. 

5.4 Pairwise Velocity Dispersion 

Pairwise velocity dispersions (PVP) are a useful tool be- 
cause of th eir strong sensitivity to the abundance of mas- 
sive haloes llMo et al.|[l99l IJing fc Bornerll2004l ; lYang et all 
12004 ITinker et alj|2006l) . Fig. [3] shows that the present-day 
abundance of massive haloes with mass > 10 14 M Q //i 111 
WMAP3 is almost 2.5 times lower than the corresponding 
abundance in WMAP1. This should be reflected in differ- 
ences in the predicted PVP. 

In Fig. O we present the dependence of PVP on 
SPSS ro.i luminosity for the three SA models used in 
this work. The black points with er ror bars show observa- 
tional estimates from lLi e t al. (2007) obtained by modelling 
the full two-dimensional redshift-space correlations. In or- 
der to compare our model prediction with these data, we 
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Figure 9. Pairwiso velocity dispersions for our three SA models 
and for four luminosity bins in the SDSS rrj.i band. Results from 
our three models are presented in blue (Al), red (B3) and green 
(C3). The black sym bols with error bars are observational data 
from (iLi et alj|2007h . The horizonal dotted line corresponds to 
the 012 = 500km/s 



adopt t he assumptions and modelling methods of lJing et al.l 
|l998t ); lLi et~aH (|2007h except that rather than construct- 
ing large mock catalogues, we measure directly both the 2- 
dimensional correlation function in redshift space (£(r p ,n)) 
and the real space correlation (£(r)) using the distant ob- 
server approximation for the former. Fig.[9]shows that model 
Al produces higher dispersions (by ~ 50 — 150km/s) on all 
scales and for all luminosity subsamples than do the two 
WMAP3 models. Except for the faintest galaxies, model Al 
also gives a b e tter fi t to the observational data. This agrees 
with lLi et all (12007) who c o mpar ed their measurements to 
the model of ICroton et al.l (|2006l ) which is very similar to 
our Al. Both WMAP3 models predict lower PVD than are 
observed, and they also seem to show a different slope. Thus, 
PVD observations se em to prefer a WMA P1 cosmology, as 
can be inferred fromljing fc Bornerl i 



Figure 10. 2-point autocorrelation functions for luminous galax- 
ies in K-band at z = 0. Solid colored curves show results for the 
three "good" models of this paper. The cyan symbols with error 
bars are predictions based on 64 small boxes (each with the vol- 
ume of the simulations of this paper) cut from th e Millennium 
Simulation catalogue of |Pe Lucia &: Blaizotl ||2007|) . Black sym- 
bols with error bars and the black solid line sho w observational 
estim ates and a power law fit from the 2DFGRS l lHawkins et al.l 
2003). The residuals shown in the lower panel are calculated with 
respect to this power law fit. 



5.5 Galaxy Clustering 

Figure \W\ compares the 2-point correlation function of our 
model galaxies at redshift z = to a recent measurement 
from the 2dFGRS j|Hawkins et al.ll2003l ). Solid colored lines 
show predictions from our three "good" models (blue for 
model Al, red and green for models B3 and C3 respec- 
tively). For each model there are about 17000 galaxies with 
Afk < —23. This magnitude limit is just fainter than the 
knee in Fig. [S] so that most luminous galaxies are included. 
In order to estimate the cosmic variance in these estimates, 
we divide the Millennium Simulation galaxy catalogue of 



157 Lucia fc Blaizotl (2007) into 64 small boxes, each with 
the same volume as our new simulations, and then calcu- 
late the mean and scatter of the resulting 64 autocorrela- 
tion estimates. (We have taken care to eliminate edge effects 
when calculating correlations for galaxies within these sub- 
volumes.) Cyan symbols show the mean calculated in this 
way, while the error bars give the 10% to 90% range. Our 
models were not tuned to match observational measures of 
galaxy clustering, so the overall agreement with observations 
demonstrated by Fig llOl is encouraging. 

We can also see that results from the three models con- 
verge on large scales (> 6/i _1 Mpc). On intermediate scales 
(6/i _1 Mpc > r > lh _1 Mpc), models B3 and C3 exhibit 
weaker clustering than model Al, agreeing better with the 
observational estimates. Interestingly, the 2-point correla- 
tion function of subhaloes (see Fig. [4| shows a similar off- 
set between the models on these scales. On smaller scales 
(r < l/i _1 Mpc) the three models give different results - 
B3 agrees with the observational data to within the scat- 
ter found among the 64 Millennium samples, while C3 is 
marginally high and Al is significantly high. Clearly cluster- 
ing on these scales is quite sensitive to details of the galaxy 
formation physics. 

It is also interesting that model Al is significantly high 
compared to the mean of the Millennium results on scales 
between 2 and 8/i -1 Mpc, even though the two simulations 
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adopt the same galaxy formation physics within the same 
cosmology. This suggests that the particular realisation of a 
L = 125/i~ 1 Mpc box used in this paper overestimates clus- 
tering on these scales. Correcting for this would bring models 
B3 and C3 into excellent agreement with the observations 
for r > 2h~ Mpc. In summary, all three models agree well 
with the data on the large scales that are sensitive to 2-halo 
correlations. Model B3 also agrees well with observation on 
smaller scales which are dominated by galaxy pairs within a 
common halo. C3 is slightly high on these scales and Al is 
too high to be compatible with the observational data. We 
note, however, that there is a significant difference between 
B3 and C3 which are implemented on the same WMAP3 
simulation. This emphasizes that small-scale galaxy correla- 
tions are very sensitive to details of the adopted galaxy for- 
mation physics and are unlikely to be useful for constraining 
cosmological parameters. 

In Fig. [11] we study how galaxy correlations vary with 
luminosity, comparing the projected autocorrelations w(r p ) 
from our models (soli d coloured line s) with observational 
data from the SDSS (|Li et al.l 120061 ') (symbols with crror 
bars). Results are shown for six magnitude bins from faint 
(top left) to bright (bottom right). Below each panel, we 
also show the ratio between the model and the observed es- 
timates. The model projected correlation function has been 
determined by integrating the real space correlation function 
(£(r)) along the line-of-sight: 




.0 10.0 
r p [Mpc/h] 



u»(r P ) = 2 / Z(Jri + r*)dr {{ =2 / g(r) 7=== (7) 

We note that if the integration is truncated at r = 60Mpc//i 
(half the box size of our simulations) the resulting projected 
correlation function is reliable up to ~ 10Mpc//i. Because 
of the limited volume of our simulations, the two bright- 
est magnitude bins contain only a few thousand and a few 
hundred galaxies respectively. Results for the fainter bins 
are based on much larger numbers of galaxies. The straight 
black line reproduced in each panel to facilitate comparison 
corresponds to the power-law: 

£ = (r/5Mpc/i _1 )" L8 

The results in this figure show reasonably good agree- 
ment between the models and the observations for M r < 
—20, but significant overpredictions, particularly for model 
Al, at fainter absolute magnitudes and at large scales. 
The differences between model and observation show sim- 
ilar trends with pair separation for all the models, and are 
as large as the differences between the models themselves. 
On the basis of this comparison none of the models is obvi- 
ously better or worse than the others. In particular, at bright 
magnitudes (M r < —20) model Al is a significantly better 
fit to the observatio ns than either B3 or C3. This reproduce 
the trends which iLi et all (|2007T l found in their own com- 
parison of S DSS correlation to t he Millennium Simulation 
catalogue of ICroton et al.l (I2006I ). This is the exact oppo- 
site of the conclusion drawn above from Fig. II 01 suggesting 
that the level of agreement between different observational 
estimates of galaxy correlations is not yet good enough to 
distinguish between the various models we are discussing. 



Figure 11. Projected correlation function w p in the SDSS 7*0.1 
band. Coloured lines show results for our si mulated galax y cata- 
logues. Black points are for the SDSS data jLi et alj|2005> . Each 
pair of panels corresponds to a different absolute magnitude bin. 
The solid black line in the upper panel of each pair corresponds to 
the power-law: £ = (r/5 Mpc/i -1 ) -18 . The lower panel of each 
pair plots the ratio of model to observation, with error bars to 
indicate the uncertainty in the observational estimate. 



5.6 Evolution to high redshift 

Figure [T2l shows how results from our models compare to the 
observed cosmic star formation rate density as a function of 
redshift. Symbols with error bars are a compilation of ob- 
servat ional data taken from Fig. 12 of ISpringel fc Hernquistl 
(2003). Models B3 and C3 have quite similar star forma- 
tion histories, although C3 lies above B3 by about 15 per 
cent at z > 2. This is due to the higher star formation ef- 
ficiency used in model C. Model Al provides much larger 
star formation rates than either B3 or C3 for z > 1 and 
lower star formation rates in the local Universe. The high 
redshift difference reflects the earlier formation of structure 
in the WMAP1 cosmology, while the difference at late times 
is a consequence of the requirement that all models produce 
the same total number of stars (as measured by the K-band 
luminosity function) in the present universe. 

The most dramatic difference in Fig.[l2]is that between 
the two cosmologies at the highest redshifts. For z ~ 6, 
the star formation rates are almost an order of magnitude 
lower in models B3 and C3 than in model Al. A related 
result, also visible in Fig. 1121 is that the peak of the cos- 
mic star formation history is shifted to lower redshift in the 
new cosmology: from z ~ 3 in model Al to z ~ 2 in model 
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Figure 12. Cosmic star formation rate density as a function of 
redshift. Symbols with error bars are a compilation of observa - 
tional data (taken from Fig. 12 of ISpringel fc Hernquistll2003l) . 
The solid coloured curves show results from our 'best' models 
(see text for details). 



B3 or to z ~ 2.5 in model C3. Fig. [12] suggests that mea- 
surements of the cosmic star formation rate at high redshift 
can potentially constrain models like those discussed here. 
Unfortunately, observational uncertainties (e.g. due to the 
use of different star formation estimators at different red- 
shifts, and to the need for substantial dust corrections) are 
too large to discriminate reliably between our three models. 
In addition, these models are far from exhausting all phys- 
ically plausible possibilities for the phenomenology of star 
formation and feedback, so the true theoretical uncertainty 
is undoubtedly larger than suggested by F ig. 1121 

In a recent paper IKitzbichler fc White] (120071) compared 



result s from the model discussed in iDe Lucia fc Blaizotl 
l|2007h to a variety of observational data at high redshift. 
They found this model (which is identical to our model A) 
to give moderately good agreement with the observed lu- 
minosity and stellar mass functions of galaxies over a large 
redshift range. As shown in Fig. 1121 the three models used 
in the present study have a significantly different behaviour 
at high redshift, so it is interesting to see if the observations 
can discriminate between them. 

Our figures 1131 and I14 | cor respond to figures 5 and 
7 of IKitzbichler fcWhitd (fcoOTl ) and show evolution with 
redshift of the rest-frame K-band luminosity function, and 
of the stellar mass function respectively. In Fig. 1131 sym- 
bols with error bars show observat i onal d eterminations 
fromlCole et~ai1 (120011 ) ; iPozzetti et all (|2003l ) ; iFeulner et all 



20031 ) ; ISaraccoet all d200ut) 



The o bserv ational estimate in 
the local Universe from lCole et al.l (|200ll ) is repeated in the 
other panels as a black dashed line. Model results are shown 
as solid coloured lines (blue for Al, red for B3 and green for 
C3). Fig. 1131 already showed all three models to agree nicely 
with local observations. At higher redshifts, the agreement 
is also reasonably good. All three seem to underpredict the 
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Figure 13. Rest-frame K-band luminosity functions for differ- 
ent redshift intervals. Model results are shown as solid lines (blue 
for Al, red for B3, and green for C3). Symbols with error bars 
show observational estimates from several surveys, as labelled 
in each panel. T he low-redshift observational determinations of 
ICole et all J200lf) is repeated as a dotted black line in the other 
panels. 



number of luminous galaxies in the highest redshift bin, but 
it should be kept in mind that the rest-frame K-band lu- 
minosities here have been extrapolated beyond the directly 
observ ed region, and so are quite uncertain (see the discus- 
sion m IKitzbichler fc White! 120071 ) . Our three models start 
to show significant differences for z > 1. Note that model Al 
lies between models B3 and C3 in this plot so that the dif- 
ferences are due mainly to galaxy formation physics rather 
than to cosmological parameters. They are, in any case, com- 
parable to the uncertainties in the observations, so that sig- 
nificantly better data are required at these redshifts to put 
strong constraints on our models. Finally, we note that the 
error bars of Fig. [T3] underestimate the true uncertainties 
as they do not include the effects of cosmic variance. These 
are particularly important when small regions of the sky are 
sampled - as is the case for the data in the highest redshift 
bin. 

Fig. [14] shows the evolution of the galaxy stellar mass 
function for our three "good" models (coloured solid and 
dashed lines) and compares them w i th observational de - 



Drory et al 



2005) 



2007), 



terminations f r om Cole et all (12001 
iFontana et all (|2006l ). As in IKitzbichler fc White! 
model results are shown both with (solid) and without 
(dashed) convolution with a normal distribution of standard 
deviation 0.25 dex, intended to represent measurement er- 
rors in logM*. The dashed black lin e in each panel rep eats 
the local observational estimate by I Cole et aLl (|200ll ). At 
redshifts beyond 1, model Al predicts substantially more 
galaxies in the mass range 9 < log(M*/M0) < 10.6 than ei- 
ther of the models in the WMAP3 cosmology. This remains 
true at higher redshift for model B3, but not for model C3. 
For log(M,/M ) > 10.8 models Al and CI predict similar 
numbers of massive galaxies at all redshifts. The increased 
star formation efficiency in model C as compared to model 
A clearly compensates for the lower as in the WMAP3 cos- 
mology. This demonstrates that even at high redshift it may 
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Figure 14. Evolution of the galaxy stellar mas s function from z = to z = 4.5. Model results are shown as solid lines (blue for Al, red 
for B3, and green for C3). Local data are from lCole et al. I ll200ll ) and are repeated as a black dashed curve in the high redshift panels. 
High redshift data are taken from Drory et al. (2005, symbols) and Fontana et al. (2006, grey shaded areas). Model predictions are shown 
both with (solid) and without (dotted) convolution with a normal distribution of standard deviation 0.25. At z = we consider the mass 
determinations precise enough to neglect this effect. 



be difficult to use galaxy data to distinguish between cos- 
mologies unless the physics of galaxy formation can be in- 
dependently constrained. Model B3 severely underpredicts 
the number of massive objects at z = 4.5, but even this dis- 
agreement may not be significant once cosmic variance and 
observational uncertainties are taken into account. 

Finally, we look at the evolution of clustering to high 
redshift. Fig. [15] presents spatial 2-point correlation func- 
tions at five different redshifts for all resolved subhaloes (for 
comparison with the lower panel of Fig. 2} and for galax- 
ies in four different stellar mass bins in each of our three 
"good" galaxy formation models. Perhaps surprisingly, al- 
though the result of section [4] that the clustering of resolved 
subhaloes is very similar in our two cosmologies, is even more 
accurately true at z = 1, it does not hold out to high red- 
shift. For z > 3, resolved subhaloes are actually substantially 
more strongly clustered in the WMAP3 cosmology than in 
WMAP1, despite the fact that the former has a significantly 
tower mass clustering amplitude. This effect is also visible in 
the galaxy autocorrelations. At the present day our models 
predict galaxies of all stellar masses to be somewhat more 
strongly clustered in the WMAP1 cosmology than in the 
WMAP3 cosmology. Beyond z — 2 the opposite is true. By 
z = 5 the effect is quite strong, more than a factor of two 
in correlation amplitude. These effects may seem surprising, 
but in fact the corresponding results for dark halos are easily 



obtained if standard analytic models are applied to our two 
cosmologies. For the convenience of the reader we provide 
an appendix repeating Mo & White's (2002) graphical anal- 
ysis of evolution in the abundance and clustering of halos for 
these two cases. The large difference in the predicted galaxy 
clustering properties in the two cosmological model could 
potentially help to diagnose cosmological parameters. We 
note, however, that our model is not able to predict reliably 
high redshift galaxy populations such as Ly-a emitters, ul- 
traluminous infrared galaxies etc. In addition, observational 
measurements at these redshifts are affected by systematics 
that are not currently well understood. A rigorous compar- 
ison between model predictions and observational measure- 
ments at these redshifts is thus still difficult. 



6 CONCLUSIONS AND DISCUSSION 

We have carried out cosmological structure formation simu- 
lations of a ACDM Universe for the cosmological parameter 
sets suggested by the first- and third-year WMAP results. 
The significant reduction in the best value for the amplitude 
of matter fluctuations on 8 ft -1 Mpc scale (as) combines with 
the decrease in the estimate of the scalar spectral index for 
primordial perturbations (n), and with the lowered matter 
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Figure 15. Autocorrelation functions for galaxies and for resolved subhaloes at five different redshifts (z = 4.9; 3.0; 2.0; 1.0; 0.). The 
galaxy results are given for four disjoint ranges in stellar mass, as indicated. (Stellar mass M* is given in units of Mq .) For M* > 3 X 10 10 
there are too few galaxies at z = 4.9 to get meaningful results, so we leave this panel blank. Different colours in the galaxy panels refer to 
our different formation models, blue for Al, red for B3, and green for C3. In the subhalo column the colours refer to our two simulations, 
blue for WMAP1, red for WMAP3. We have multiplied all correlation functions by r 18 to make the differences between the models more 
visible. Note that r is a comoving coordinate with units of h~ 1 Mpc. 



density (fl m ) to produce a significant delay in structure for- 
mation in the WMAP3 case (see Sec. QJ . 

By coupling our numerical simulations to semi-analytic 
models for galaxy formation, we have investigated the im- 
plications of this delay for the observed properties of galax- 
ies, both at low and at high redshift. Specifically, we 
have compared the galaxy formation model described in 



|Pe Lucia fc BlaizotJ (|2007t) for the WMAP1 cosmology to 
two galaxy formation models for the WMAP3 cosmology 
which use the same physical framework but different effi- 
ciency parameters. We find that both new parameter sets 
can compensate for the delay in structure formation to pro- 
duce galaxy populations at z = which agree with observa- 
tion just as well as the old model for the WMAP1 cosmology. 
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The luminosity functions are almost identical, the correla- 
tion functions show at most small differences, and offsets in 
the predicted Tully-Fisher relations are difficult to interpret 
because disk rotation velocities cannot be predicted reliably 
to the level of accuracy required. 

Pairwise velocity dispersion measurements are sensitive 
to cluster abundance and therefore differ significantly for the 
two cosmological models of our study. We have shown that 
the two WMAP3 models underpredict the measured PVD 
by more than fOOkm/s on scales 0.3 < r p < 2Mpc/h, and 
are lower than the corresponding predictions of the WMAPf 
model by 50 ~ f50km/s on all scales. 

Substantial differences between the various models ap- 
pear at high redshift. The delay in the structure formation 
translates directly into a delay in the global star formation 
history: at z ~ 6 the star formation rates in the models 
based on WMAP3 are lower than those based on WMAPf 
by almost an order of magnitude. As discussed elsewhere, 
this has important implications for the formation of the first 
stars, and for reionization. Predictions of our three mod- 
els for galaxy luminosity and mass functions at high red- 
shift show substantial differences. Unfortunately, the uncer- 
tainties in the observed luminosities and masses, combined 
with large cosmic variance uncertainties, are still too large 
to place strong constraints on the efficiencies and scalings 
of the physical processes we model. Somewhat counterintu- 
itively, we find that at high redshift galaxies of given stellar 
mass are predicted to be substantially more clustered in the 
WMAP3 cosmology than for WMAPf. 

When comparing z — correlation functions from our 
three models to recent observatio nal determinati ons from 
the 2 dFGRS (|Hawkins et aUl2003h and the SDSS (|Li et all 
2006), we found interesting and apparently significant differ- 
ences between the models, particularly for galaxies fainter 
than ~ —20 in the SDSS r-band. At these magnitudes, our 
WMAP3 galaxy catalogues provide better agreement with 
the SDSS data but the opposite is true for galaxies around 
the knee of the luminosity function. It is important to re- 
alise that the reduced mass clustering amplitude implied by 
the WMAP3 parameters is almost entirely offset by an in- 
crease in halo bias, so that predictions for galaxy clustering 
change very little. At least at separations r < 20/i _1 Mpc, 
galaxy clustering is much more sensitive to galaxy forma- 
tion physics than to cosmological para meters (see also the 
discussion in Ivan den Bosch et al J 12003). In fact, for almost 
all of the population properties we have examined, the vari- 
ations induced by "acceptable" variations in the galaxy for- 
mation parameters are at least as large as those produced by 
the variation in cosmological parameters between WMAP1 
and WMAP3. The systematic properties of galaxies and 
their small-scale clustering should be used to understand 
how galaxies form, not to constrain cosmolo g y. A s imilar 
conclusion was reached by iKauffmann et alJ (jl999ah who 
used cruder semianalytic models to show that similar z = 
clustering was predicted in two quite different cosmologi- 
cal models (rCDM and ACDM). In this case, however, the 
predicted evolution to high redshift was d ifferent enough to 
offer a clear way to distinguish the models (|Kauffmann et ah! 
Il999bh . 

In this paper we have shown that varying efficiency pa- 
rameters within a given framework for modelling galaxy for- 
mation can lead to very similar predictions for the evolution 



and clustering of galaxies in the WMAP1 and WMAP3 cos- 
mologies. The problem is sufficiently degenerate that a va- 
riety of acceptable parameter sets can be found in either 
cosmology. It may seem unsatisfactory to 'fine-tune' model 
parameters to fit the observational data, but it is interesting 
that substantially different efficiencies of star formation and 
feedback are required in the various cases. More detailed ob- 
servational data on how these processes work in individual 
systems may therefore shed light on which parameter values 
are appropriate. In addition, differing efficiencies translate 
into significantly different predictions at high redshift. More 
detailed and more statistically precise observations of high- 
redshift galaxies will be able to distinguish between the mod- 
els. Once these aspects of the galaxy formation process are 
better understood, it may indeed be possible to use galaxy 
surveys to constrain cosmological parameters. 
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Figure Al. The evolution of the comoving number density of 
dark matter haloes with mass exceeding a specific value M in our 
two cosmological models: WMAP1 (black curves) and WMAP3 
(red curves) . The numbers labelling the black curves indicate the 
value of log(M/MQ) for the WMAP1 case. See the text for more 
details. 
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APPENDIX A: HALO ABUNDANCE AND 
CLUSTERING 

In this appendix we use the formulae and the graphical pre- 
sentation scheme of iMo fc White! (|2002h to illustrate how 
the evolution of halo abundance and clustering differs be- 
tween our WMAP1 and WMAP3 cosmologies. As in figures 
1 and 2 of IMo fc White! (I2002T ). we present plots of comov- 
ing abundance against redshift for halo samples defined by 
lower mass limits Mmm ( z ) corresponding to a variety of halo 
properties, in particular, for lower limits which correspond 
at all redshifts to a given halo mass, a given halo virial tem- 
perature, a given fraction of the total cosmic mass density, 
a given clusering strength in comoving units, and a given 
clustering stgr ength at z = for the halo descendents. We 
refer to IMo fc White! (2002) for detailed discussion of these 

quantities and for the relevant formulae. 

Figure IA1I corresponds to figure 1 of IMo fc White! 
(2002). Each curve gives the comoving abundance as a func- 
tion of redshift of halos more massive than the value indi- 
cated by the label. Curves for our WMAP1 parameters are 
shown in black while curves for our WMAP3 parameters 
are shown in red. Labels give the decimal logarithm of halo 
mass in units of solar masses and are placed next to the cor- 
responding WMAP1 curve. Curves are plotted for the same 
set of mass limits in the two cosmologies, and give almost 
identical abundances at z = for low masses. This can be 
used to identify the WMAP3 curves at high redshift where 
they deviate very substantially from the WMAP1 curves. 



ISpringel et ail (|2005l ) showed these theoretical predictions 
to be in excellent agreement with Millennium Simulation 
results for z ^ 10. At all masses the difference in halo abun- 
dance between the two cosmologies increases with the in- 
creasing redshift. For M > 10 M© the abundance differ- 
ence is already almost an order of magnitude at z = 0, and 
the same is true for M > 10 12 M Q at z = 5 and M > 1O 8 M 
at z = 10. 

In Fig. lA2l we give abundance-redshift relations for halo 
samples defined above lower mass limits M m i n (z) which cor- 
respon d to other conditions, as in figure 2 of IMo fc White] 
(2002). In each of these plots the two curves of Fig. I All are 
repeated as dashed curves. These can be used as a reference 
to obtain the halo mass corresponding to each point in the 
abundance-redshift plane. 

In the upper left panel of Fig. IA2I the solid curves link 
halo populations containing given fractions F of the total 
cosmic mass density at each redshift. The labels give F val- 
ues for the curves they are placed next to. Where no label 
is given the F value can be inferred from the surrounding 
curves. Red curves show WMAP3 results for the same F 
values and lie above the corresponding WMAP1 curves at 
all redshifts. At z = the curves are close enough that it 
is easy to infer the F value for each red curve by compar- 
ing it with the corresponding black curve. For example, for 
WMAP1 at z = 0, one percent of all cosmic mass is in dark 
halos above a lower mass limit corresponding to abundance 
n = 5 x 10~ 7 /i 3 Mpc" 3 , thus M > 1O 15 M (from Fig. |A"TJ) . 
For WMAP3, the z = abundance at F = 0.01 is about 
a factor of 2 higher, and the corresponding mass is about 
2.5 times smaller. At z = 5 the one percent mass point cor- 
responds to n = 10~ 3 /i 3 Mpc~ 3 and M > 4 x IO^Mq for 
WMAP1, but to n = 10" 2 /i 3 Mpc" 3 and M > 4 x 1O 1O M 
for WMAP3. 

In the upper right panel of Fig. IA2l the solid curves link 
halo populations at each redshift with virial temperatures T 
in excess of a given value. Labels give the decimal logarithm 
of the limiting temperature in Kelvin and are placed next to 
the curve they refer to. At z — 0, there is a close correspon- 
dance between WMAP1 and WMAP3 at high abundance. 
At low abundance (high mass) the WMAP3 curves lie be- 
low their WMAP1 counterparts. At z = 10, the current 
best estimate of the reionization redshift, only halos with 
M > 8 x 1O 7 M have virial temperatures sufficient to ion- 
ize hydrogen (T > 10 4 K) and are thus able to cool their 
baryonic component effectively. For WMAP1 parameters, 
the comoving abundance of such halos is n = 10/i 3 Mpc" 3 
and they contain a fraction F ~ 0.04 of all cosmic matter 
(from the upper left plot of Fig. lA"2")l . For WMAP3 parame- 
ters, the predicted abundance of such halos drops by about 
a factor of 6, and the fraction of cosmic matter contained 
in them drops by about an order of magnitude. This is the 
reason why reionization is much more difficult to explain (at 
z — 10) for the revised WMAP parameters. Differences of 
this kind also explain why we get higher global star forma- 
tion rates at z ~ 5 in our model Al than in our models B3 
and C3 (Fig|12p as well as a correspondingly higher stellar 
mass function in Al at high redshift (Fig. I14|l . 

In the lower left panel of Fig. IA2I the solid curves link 
halo populations at each redshift which have a given strength 
of clustering in comoving coordinates as characterized by 
Ag(M, z), the rms fluctuation in overdensity of haloes more 
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Figure A2. A repeat of figure 2 of lMo fc White! ||2002T ) but showing results for both WMAP1 (black) and WMAP3 (red) parameters. 
The limiting mass of the halo population is chosen so as to keep a different quantity constant along a solid curve in each panel: cosmic 
mass fraction F (top left), minimum virial temperature T (top right), clustering strength Ag (bottom left) and clustering strength of 
the z = descendents As,o (bottom right). The dashed curves in each panel repeat those of Fig. lAll The black and red numbers label 
the WMAP1 and WMAP3 curves, respectively. See the text for details. 



massive than M at redshift z after smoothing with a spher- 
ical top-hat filter of comoving radius 8/i _1 Mpc. As in the 
other panels, the labels give As values for the curves they 
are placed next to. At redshifts 5 to 10, the WMAP3 curves 
lie below the corresponding WMAP1 curves at low abun- 
dances but above them at high abundances. Comparing the 
dashed and solid curves in this panel, one sees that at z — 
the clustering strength of halos is very similar in the two 
cosmologies for all lower limits to halo mass. This is the re- 
sult visible in the lower panel of Fig. [4] At z = 5, however, 
these same curves show the clustering strength of halos to 
be substantially stronger in WMAP3 for all limiting halo 
masses than in WMAP1 (typically by 25 to 40% in A 8 ). 



This is the result seen for low-mass (sub)halos in the right 
column of Fig. 1151 

Finally, in the lower right panel of Fig. IA2l solid curves 
link halo populations at each redshift for which the z — 
descendents have a given clustering strength, as indicated 
by Ag,o, their present-day value of As. As in the other pan- 
els, the labels give values of As,o for the curves they are 
placed next to. At z = 0, the two cosmologies predict nearly 
identical clustering strengths at intermediate abundances. 
iMo fc White! (|2002r ) give a number of examples of how this 
plot may be used. Here we note that the progenitor popula- 
tion which ends up with a given z = clustering strength is 



© 2007 RAS, MNRAS 000,ITlfl"4l 



18 J. Wang et al. 



substantially more strongly clustered at high redshift in the 
WMAP3 cosmology than in the WMAP1 cosmology. 

This paper has been typeset from a TfjjX/ KTfijX file prepared 
by the author. 
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